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^^ ■ Abstract. Prior to the epoch of full reionization, the intergalactic medium 

and gravitationally collapsed systems will be detectable in 21-cm radia- 
r^ • tion. Physical mechanisms that would produce a 21-cm signature are 

[t . discussed. These include Lya coupling of the hydrogen spin temperature 

to the kinetic temperature of the gas resulting from the radiation by an 
early generation of stars, preheating by soft x-rays from collapsing dark 
matter halos, and preheating by ambient Lya photons. A patchwork of 
,-H ■ either 21-cm emission, or absorption against the Cosmic Microwave Back- 

^ I ground, will result. Large radio telescopes like the Giant Metre Wave- 

''^ . length Telescope (GMRT) or a Square Kilometre Array Radio Telescope 

_^ I (SKA) offer the prospect of measuring this signature, and so detecting 

^\i . the transitional epoch from a dark universe to one with light. 

o: 
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^ ■ 1. Introduction 

C^ ■ 

K^ ■ The development of structure in the Universe was well advanced at early times. 

k^ \ Quasars have been detected nearly to a redshift of z = 5, and the most distant 

\^ • galaxies to even greater redshifts (Dey et al. 1998). The spectra of high redshift 

QSOs have additionally shown that the Intergalactic Medium (IGM) itself had 
undergone an extensive development of nonlinear structures at early times as 
well, as revealed by the Lya forest. Still unknown, however, is the nature of 
energetic processes at these early times. While numerical simulations have shown 
that the IGM is expected to fragment into structures at early times in Cold Dark 
Matter (CDM) dominated cosmologies (Zhang et al. 1998), and even into early 
galaxies (Governato et al. 1997), the simulations are much less able to predict 
the efficiency with which gravitationally collapsed objects will emit radiation. 
Although QSO sources may account for the photoionizing UV background at 
high redshifts (Meiksin & Madau 1993; Haardt & Madau 1997), it is less clear 
that they were responsible for the original reionization of the IGM. Similarly, 
although z > b galaxies have been detected, the epoch during which the first 
generation of stars formed is still poorly constrained (Madau et al. 1998; Hughes 
et al. 1998). Although IR observations will permit even higher redshift galaxies 
and QSOs to be observed, detections become increasingly difficult because of the 
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diminution in surface brightness due to cosmological expansion. This difficulty 
cahs for alternative means to be found for discovering the nature and period of 
the first major generation of energy-producing sources, the epoch of First Light. 
I discuss here possible ways in which a Square Kilometre Array (SKA) radio 
telescope could reveal the first epochs of energy generation. Because so little is 
known about the distribution and nature of such sources, the emphasis will be 
on the physical mechanisms related to their detection. The means by which the 
sources are revealed is through their impact on the surrounding neutral IGM 
and the resulting emission or absorption of 21-cm radiation. In the next section, 
I review the Wouthuysen-Field mechanism for the excitation of 21-cm emission 
through the scattering of Lya photons, and the heating of the neutral IGM 
by resonant scattering of the Lya photons and soft x-ray heating from early 
collapsed objects. In §111, a few scenarios are presented for the possible effects 
of early generations of radiation sources on the IGM and their detection. 

2. 21-cm Emission and Absorption Mechanisms 

2.1. The Spin Temperature 

The emission or absorption of 21-cm radiation from a neutral IGM is governed 
by the spin temperature Ts of the hydrogen, defined by 

— = 3exp[-r*/Ts], (1) 

no 

where no and rii are the singlet and triplet n = 1 hyperfine levels, T* = 
hi^io/kB = 0.07 K, where viq is the frequency of the 21-cm transition, h is 
Planck's constant, and k^ is Boltzmann's constant. In the presence of only the 
Cosmic Microwave Background (CMB) radiation, the spin temperature will be 
the same as the temperature of the CMB, and no emission or absorption relative 
to the CMB will be detectable. A mechanism is required that decouples the two 
temperatures. This may be achieved by coupling the spin temperature to the 
kinetic temperature of the gas itself. Two mechanisms are available, collisions 
between hydrogen atoms (Purcell & Field 1956) and scattering by Lya photons 
(Wouthuysen 1952; Field 1958). The collision-induced coupling between the 
spin and kinetic temperatures is dominated by the spin-exchange process be- 
tween the colliding hydrogen atoms. The rate, however, is too small for realistic 
IGM densities at the redshifts of interest, although collisions may be important 
in dense regions, 6p/p ;^ 30[(1 + z)/lO]-'^ (Madau, Meiksin, & Rees 1997). 

Instead the dominant mechanism is likely to be Lya scattering through 
the Wouthuysen-Field effect. This process mixes the hyperfine levels of neutral 
hydrogen in its ground state via an intermediate transition to the 2p state. An 
atom initially in the n = 1 singlet state may absorb a Lya photon that puts 
it in an n = 2 state, allowing it to return to the triplet n = 1 state by a 
spontaneous decay. At this point, the astute student of quantum mechanics 
will ask how is it possible for electric dipole radiation (Lya photons) to induce a 
spin transition? The key is spin-orbit coupling: it's the total angular momentum 
F = I + J that counts. (Here / is the proton spin and J is the total electron 
angular momentum, J = S + L.) There are four hyperfine states involved, the 
n = 1 singlet O'S'1/2 Siiid triplet i^i^ states (the notation is fLj), and the two 
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triplet n = 2 states 1P1/2 and 1-P3/2- The selection rule AF = 0, 1 permits the 
transitions 0S1/2 -^ 1P1/2, 1P3/2 and 1P1/2, 1-P3/2 ^ 15*1/2, and so effectively 
oS'1/2 -^ i5'i/2 occurs via one of the n = 2 states. 

When the IGM is highly opaque to the scattering of Lya photons, as it is 
when still neutral, the large number of scatterings of Lya photons in an ambient 
radiation field will ensure a Boltzmann distribution for the photon energies near 
the Lya frequency, with a temperature given by the kinetic temperature Tk of 
the IGM (Field 1959). In this case, the spin temperature of the neutral hydrogen 
becomes^ 

^ TcMB + UaTK /„x 

1 + ya 

where Tcmb = 2.73(1 + z) K is the temperature of the CMB (Mather et al. 
1994), and 

Va = -^—7^ (3) 

is the Lya pumping efficiency. Here, ^10 = 2.9 x 10~^^s~^ is the spontaneous 
decay rate of the hyperfine transition of atomic hydrogen, Pio is the indirect 
de-excitation rate of the triplet via absorption of a Lya photon to the n = 2 
level, and Ts ^ T* was assumed. In the absence of Lya pumping the spin 
temperature goes to equilibrium with the 21-cm background radiation field on 
a timescale T*/(rcMB^io) ~ 5 x 10^ yr, and neutral intergalactic hydrogen will 
produce neither an absorption nor emission signature. If ya is large, Ts -^ Tx, 
signifying equilibrium with the matter. A consideration of the net transition 
rates between the various hyperfine n = 1 and n = 2 levels above shows that 
the 1^0 transition rate via Lya scattering is related to the total rate Pa by 
Pio = 4Pq,/27 (Field 1958). This relation and equation (|2|) are derived in the 
Appendix. In the limit Tk ^ Tcmb, the fractional deviation in a steady state 
of the spin temperature from the temperature of the CMB is 



Ts — Tqmb 
T's 



1 + 



2cMB 



VoTk 



-1 



(4) 



There exists then a critical value of P^ which, if greatly exceeded, would drive 
Ts -^ Tk- This thermalization rate is (Madau et al. 1997) 



_ 27^ioTcMB 

-fth = 



(5.3 X 10-12 s-i) (i±i) . (5) 



4T, 

2.2. 21-cm Emission Efficiency 

To illustrate the basic principle of the proposed observations, consider a region 
of neutral material with spin temperature Ts 7^ Tqmb, having angular size on 
the sky which is large compared to a beamwidth, and radial velocity extent due 



^In the presence of a radio source, the antenna temperature of the radio emission should be 
added to Tcmb in equation 0. The radio emission may make an important contribution in the 
vicinity of a radio-loud quasar (Bahcall & Ekers 1969), and would itself permit the IGM to be 
detected in 21-cm radiation. 



to the Hubble expansion which is larger than the bandwidth. Its intergalactic 
optical depth at 21(1 + z) cm along the line of sight, 
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_ 3c3/l3nHl(0)v4io 15 _ in-2-9?,-l /^ ^CMb \ / ^IGM/^io \ . ^ 

(6) 
will typically be much less than unity. The experiment envisaged consists of 
two measurements, separated in either angle or frequency, such that one mea- 
surement, the fiducial, detects no line feature, either because there is no H i or 
because Ts '^ TcmBj and the second at Ts 7^ Tcmb- Since the brightness temper- 
ature through the IGM is T^ = Tcmbc"^ + Ts{1 — e~'^), the differential antenna 
temperature observed at the Earth between this region and the CMB will be 



6T = {l + z)-\Ts-TcMB){l-e-n - (0.011 K)h-,^ (^^f^) (^ 
where the 21-cm radiation efficiency is defined as 



1/2 
(7) 



^ = ^™ \ — Ys — j • ^^^ 

Here xm refers to the neutral fraction of the hydrogen in the region for which 
Ts / TcMB- As long as Ts is much larger than Tcmb (hence if there has been 
significant preheating of the intergalactic gas), r/ -^ xhi, and the IGM can be 
observed in emission at a level which is independent of the exact value of Ts- 
By contrast, when Tcmb ^ Ts (negligible preheating), the differential antenna 
temperature appears, in absorption, a factor ~ Tqub/Ts larger than in emission, 
and it becomes relatively easier to detect intergalactic neutral hydrogen (Scott 
& Rees 1990). 

2.3. Preheating the IGM 

The role of the spin temperature is manifest in eq. (|8|): when Ts < Tcmb the IGM 
absorbs 21-cm radiation from the CMB, while for Ts > Tcmb the IGM emits 21- 
cm radiation in excess of the CMB. In the absence of decoupling mechanisms, 
Ts = Tcmb- The presence of Lya photons with sufficient intensity will thus 
enable the IGM to be "seen." The adiabatic expansion of the Universe will 
generally bring the kinetic temperature of the IGM well below the temperature of 
the CMB. Coupling Ts to Tk will permit the IGM to be detectable in absorption. 
If there are sources of radiation that heat the IGM, however, it may be possible 
instead to detect the IGM in emission. 

Possible heating sources are soft x-rays from an early generation of QSOs or 
thermal bremsstrahlung emission produced by the ionized gas in the collapsed 
halos of young galaxies. In CDM-dominated cosmologies, the latter may be 
in sufficient number to heat the IGM above the CMB temperature by z ~ 7 
(Madau et al. 1997). 

While photons just shortward of the photoelectric edge are absorbed at the 
ionization front generated by a QSO source, photons of much shorter wavelength 
will be able to propagate to much greater distances. Most of the photoelectric 
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heating of the IGM by a QSO is accomphshed by soft x-rays. The time required 
for the radiation at the hght front to heat the intergalactic gas to a temperature 
above that of the CMB is typically 10% of the Hubble time. The H ii region 
produced by a QSO will therefore be preceded by a warming front. Note that, 
as the X-ray-heated bubbles around QSOs will survive as fossils even after the 
quasar has died, several generations Ng of quasars may actually be responsible 

-\ In 

for preheating the entire IGM. For a typical QSO age of tg ~ 3x lO^A'g yr, the 
required QSO comoving space density to heat the entire IGM to a temperature 
above that of the CMB by z ~ 6 is ~ 10"^'' Mpc~ A''"^. By comparison, the 
comoving space density of bright QSOs at z = 4 is ~ lOOA'g times larger (Warren, 
Hewett, & Osmer 1994). If all bright galaxies undergo a quasar phase, QSOs 
must have a very short lifetime, and Ng ~ 100. Soft X-rays from a few bright 
QSO sources could then prevent collapsing structures, such as protoclusters 
while still in the linear regime, from being detected in 21-cm absorption against 
the CMB. 

An additional heating source is the Lya photon scattering itself. The aver- 
age relative change in a Lya photon's energy E after having been scattered by 
a hydrogen atom at rest is 

,/S.E, hVn o , , 

{-Et) = ^ « -10-^ 9 

E muC^ 

where mn is the mass of the hydrogen atom. (It should be noted that this is 
an approximation valid only for hva ^ kT^. In the opposite limit, energy will 
flow from the atoms to the photons.) Through recoil, energy is transferred from 
photons to atoms at a rate 

A E 

Ea = -{^)hUaPa. (10) 

where Pa is the Lya scattering rate per H atom. In the case of excitation at the 
thermalization rate Pn^, equation (ITOl) becomes 



E,^ = il^.!^^li2^^(220KGyr-^){^). (11) 



4 TTT-HC^ 1* \ 7 

(Madau et al. 1997). The characteristic timescale for heating the medium above 
the CMB temperature via Lya resonant scattering at this rate is 

about 20% of the Hubble time at z ~ 8. The result is a finite interval of time dur- 
ing which Lya photons couple the spin temperature to the kinetic temperature 
of the IGM before heating the IGM above the CMB temperature. If Lya sources 
turned on at redshifts z^ <^ 10, this interval would present a window in redshift 
space near z ~ 8 that would enable a large fraction of intergalactic gas to be 
observable at ~ 160 MHz in absorption against the CMB, and so isolate the 
epoch of First Light. 



3. Scenarios 

The calculations in this section were done in collaboration with Paolo Tozzi and 
Piero Madau. 

Two cosmological models are considered, a tilted Cold Dark Matter model 
(tCDM) with Qq = 1, Hq = 50kms^^Mpc~^, and as = 0.55, designed to match 
both CMB measurements on large scales and the constraint on amplitude im- 
posed by galaxy cluster abundances on small, and a flat open CDM model 
(OCDM) with 1^0 = 0.3, J^a = 0.7, Hq = 70kms~iMpc~\ and ag = 1.1, which 
similarly matches both constraints. The baryon density in both models is as- 
sumed to be ^Bh"^ = 0.024, where h = i7o/100kms"^Mpc~^ 

We consider a scenario in which sources of Lya photons are in sufficient 
abundance throughout the universe to couple the spin temperature to the kinetic 
temperature of the IGM everywhere. We further suppose the IGM has been 
preheated to a temperature well above that of the CMB, either by the same 
Lya photons responsible for the coupling or by soft x-ray sources. In this case, 
the IGM will emit at a rate independent of Ts (the hyperfine levels will be 
occupied according to their statistical weights, ni/riQ = 3). Because of structure 
in the IGM, the emission will not be uniform. 

In Figure 1, we show the range of density fluctuations that would be de- 
tectable in a single beam by a Square Kilometre Array Interferometer (SKAI) at 
160 MHz {z = 8), as a function of beam size and frequency band width, and for 
several assumed integration times. The detection thresholds are scaled according 
to a continuum rms in a 80 MHz band at 160 MHz of 64 nJy over an 8 hour in- 
tegration time (Braun 1998; also see www.nfra.nl/skai/science/index.htm). The 
experiment is based on taking differences between beams. The dashed lines 
show curves of constant rms antenna temperature fluctuation within the IGM, 
(((5Ta)^)^'^ = (JpTa, where ap is the rms relative density fluctuation of the IGM 
for a volume corresponding to a given bandwidth (Az//z/ = Az/[1 + z]) and 
angular size A9. Here, T^ is the mean antenna temperature from the IGM at 
redshifted 21-cm. At a fixed bandwidth, the antenna temperature fluctuation 
increases with decreasing angular scale because ap increases with decreasing lin- 
ear scale. Because the detected flux is proportional to the solid angle of the 
beam, the detected signal decreases with decreasing angular scale until it falls 
below the detection threshold, indicated by the solid lines. The size of the de- 
tectable IGM fluctuations differs greatly for the two models. This is because the 
growth of density fluctuations ceases early on in an open universe, so that the 
fluctuations on a given angular scale are much larger in the OCDM model than 
in tCDM at high redshift. 

In a second scenario, the spin temperature is again coupled to the IGM 
kinetic temperature everywhere, but the IGM has not had time to heat above 
the CMB temperature. We then consider the emission signature resulting from a 
QSO soon after it turns on as the medium surrounding it is heated by soft x-rays 
from the QSO. The experiment in this case is done by differencing beams pointed 
in regions around the QSO with one pointed through the QSO H ll region, where 
both emission and absorption are absent. An image of the resulting emission for 
the tCDM model is shown in Figure 2. The model was computed using Hydra 
(Couchman, Thomas, &: Pearce 1995). As the warming front produced by the 
QSO expands, a growing amount of the surrounding IGM is revealed. Note that, 
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Figure 1. Detectability of fluctuations in 21-cni brightness tempera- 
ture (10^2 < {{STaf)'^^'^ < 10^2), from the IGM by SKAI as a function 
of beam width and bandwidth. The dashed hues show the size of the 
rms fluctuations expected in a volume defined by the angle and fre- 
quency widths. The solid lines show the threshold below which the 
fluctuations are undetectable for integration times of 10, 100, and 500 
hours. 



although the QSO was placed in the corner of the simulation volume, the figure 
can equally be viewed as the emission due to heating by a beam of soft x-rays 
from the QSO with an opening angle of 90°. Thus imaging the gas surrounding a 
QSO in 21-cm emission would provide a direct means of measuring the opening 
angle of QSO emission. 

It should also be noted that a region of a given fixed density fluctuation 
will not always yield the same fluctuation in 21-cm emission. This is because 
of the dependence of the spin temperature on the temperature of the IGM. 
Only when the IGM temperature much exceeds that of the CMB will the 21-cm 
emission be independent of the IGM temperature, according to equation (g). 
In general, the fluctuations in brightness temperature will depend on both the 
density fluctuations of the IGM and the temperature fluctuations, which in 
turn depend on the ages and distribution of the sources. A knowledge of the 
cosmological density fluctuation spectrum, as may be measured by future CMB 
missions like MAP and Planck, will then enable the statistical distribution of the 
sources that heat the IGM, whether QSOs as here or Lya photons from early 
stars, to be established using measurements of the fluctuations in the 21-cm sky 
at high redshift. 
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A The Wouthuysen-Field Mechanism 

There are six hyperfine states in total that contribute to the Lya transition. 
These are shown in Figure 3. Only four of these, the two n = \ states and 
the two n = 2, F = 1 states, contribute to the excitation of the 21-cni line by 
Lya scattering. Denoting the occupation number of o5'i/2 by uq and that of 
i'S'i/2 by ni, the rate equation for uq is 

^ = Aio (l + f^) ^1 - 3^10 J^o + Pfo^i - Poino, (13) 

where the radiation intensity at the 21-cm frequency has been expressed in 
terms of the antenna temperature Tr, Iiq = 2k^TYi/X\Q. The ratio T^/T^ is 
the number of 21-cm photons per mode. Here, Pqi and P^q are the effective 
excitation and de-excitation rates due to Lya scattering. These may be related 
to the total scattering rate of Lya photons. Pa, as is now shown. 

It is convenient to label the n = 2 levels as states 2 — 5, from lowest energy to 
highest. Then the external radiation intensity at the frequency z/jj corresponding 

to the i ^> j transition may be expressed as the antenna temperature Tp', as 
above. The temperature corresponding to the energy difference itself, hvij, is 
expressed as Tij = hvij/k-Q- If Aji denotes the spontaneous decay rate for the 
transition j ^^ i, then Pq^ and P{q are given by 

To3 A^Q + A31 ro4 A4Q + A41 



and 



pa ^ ^13 ^30^31 T^i ^40^41 .. ^^ 

10 Ti3 .430 + A31 Tu ^40 + ^41 ■ ^ ^ 



p 

2 3/2 



1 1/2 



0^ 1/2 




p 

13/2 



p 

11/2 



p 

1/2 



Figure 3. The hyperfine energy level diagram for the Lya transition. 
The two n = 1 states are on the left, and the n = 2 states are on the 
right. The dashed lines represent forbidden transitions (F = — > F = 
and AF = 2 are not permitted), while the dot-dashed lines represent 
complementary transitions that give no contribution to the excitation 
or de-excitation of the 21-cm transition. 
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The ratios Aji/A^, where A^ is the total Lya spontaneous decay rate, may 
be solved for using a sum rule for the transitions. This states that the sum of all 
transitions of given nFJ to all the n'J' levels (summed over F' and mpi) for a 
given n'J' is proportional to 2F+1J2I There are four sets of downward transitions 
to n' = 1, J' = 1/2, corresponding to the decay intensities /51, /5o(= 0), /41, 
-^40, ^31, ^30, hi, and /2o(= 0). These give 

-^51 5 ho + hi . -^30 + -^31 o rir\ 

- - 1, = 3. (16) 



/40 + hi 3 I30 + /31 hi 

Similarly, there are four sets of upward transitions to n' = 2, J' = 1/2 or 3/2, 
giving the intensities /30, hoi= 0), /31, hi, hoi= 0), /40, -^51, and hi- These 
ratios are 

-^40 ^ 1 -^30 ^ 1 Q„x 

hi + -^51 3 hi + -^31 3 

Using hj/h = {gk/gtot)iAkj/Aa), where h is the total Lya decay intensity 
summed over all the hyperfine transitions, and g^ = 2F + 1 is the statistical 
weight of level k {gtot is the sum of the weights of the upper levels), gives 
A2o/Aa = A^q/A^ = 0, A2i/Aa = A^xjA^ = 1, A^q/A^ = A^/A^ = 1/3, 
and A3i/Aa = A^o/Aa = 2/3. We then obtain P^ = 3Aa{TR/Ta) (T„ = 
hua/kB), and Pio = (4/27)Pq,. Statistical equilibrium for Tr = requires 
P01/-P10 = ni/no = 3exp(— T*/Tft:). Solving the rate equation ([l^ for Tr > 
gives equation (§), where Tr = TcmBi and where Ts » h, Tk >> T*, and 
Tr >> r^, have been assumed. 



See, e.g., Bethe, H. A., Salpeter, E. E. 1957, Quantum Mechanics of One- and Two-Electron 
Atoms, New York: Academic Press. 

11 



